Evolution of the Spiral Waves in Excitable System 
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Spiral wave, whose rotation center can be regarded as a point defect, widely exists in various two 
dimensional excitable systems. In this paper, by making use of Duan's topological current theory, 
we obtain the charge density of spiral waves and the topological inner structure of its topological 
charge. The evolution of spiral wave is also studied from the topological properties of a two- 
dimensional vector field. The spiral waves are found generating or annihilating at the limit points 
and encountering, splitting, or merging at the bifurcation points of the two-dimensional vector field. 
Some applications of our theory are also discussed. 

PACS numbers: 03. 65. Vf, 02. 10. Kn, 82. 40. Ck 



I. INTRODUCTION 

Rotating of spiral waves is a fascinating example of self- 
organization in various two dimensional excitable sys- 
tems including Belousov-Zhabotinskii reaction [J, 0, S| 
, cardiac muscle [3, and CO oxidation on platinum 
surfaces^. Recently, much effort has been put into the 
development of the analytic description of excitable sys- 
tem, and the spiral waves have drawn great interest and 
have been studied intensively in many ways [1, 0, H, H, 

ini. 

The rotation centers of the spiral wave can be regarded 
as a point defect and defined in terms of a phase singular- 
ity. The control of these singular points is a very complex 
problem in all two-dimensional excitable system. The dy- 
namics of spiral wave are determined not only by prop- 
erties of the excitable system but also by the topological 
characteristic number of the spiral wave. The topological 
charge, which describes the strength of spiral wave phase 



plays an essential role during the processes of creating 
or annihilating new spiral waves. These processes can 
be regarded as topological phenomena in excitable sys- 
tem. The mathematical description of these spiral waves 
relates to the notion of phase which in turn allows one 
to characterize spiral waves by an index. The phase field 
9{r, t) of spiral waves is a continuously differentiable func- 
tion except at the singular points 03, |^,[2^. The sites 
of spiral wave cores are just the phase singular points. An 
index, which counts the winding number of phase field 
around the singular point, is just the topological charge 
of the spiral wave. From this description, a number of 
topological arguments have been applied to help us to 
understand the topological properties of spiral wave and 
some important topological constrains on behaviors of 
spiral wave have been investigated^, [H, [H, [H, 



ll8|, lia, 120, |2l|, iM However, the tradi- 

tional topological theory of spiral wave is not so powerful 
at directly computing the singularities of phase fields, but 
only considers a closed path around the studied region, 
and uses the single-value condition to indirectly analyze 
the existence of spiral wave in the region. Therefore, it 
is very important to develop a new method which can 
directly deduce these spiral waves in a precise topolog- 
ical form, and can reveal the inner structures and the 
branch processes (generate, annihilate, encounter, split, 
or merge) of the spiral waves. 

In this paper, by making use of Duan's topological cur- 
rent t/ieor?/[23, da, |29, 30, 3l|, we study the inner struc- 
ture of spiral wave in details. It is showed that the center 
of spiral wave is not only classified by topological charge, 
but also classified by Hopf index, Brouwer degree in more 
detail in geometry. We also investigate the evolution of 
spiral wave. And one sees that the spiral waves generate 
or annihilate at the limit points and encounter, split, or 
merge at the bifurcation points of the spiral wave field. 
Furthermore, we summarize the topological constrains on 
the behaviors of spiral waves during the above processes. 
Some possible applications of our topological theory of 
spiral wave are also discussed. 



II. THE TOPOLOGICAL CHARGE DENSITIES 
OF SPIRAL WAVES 

Excitable system is typically described in terms of 
reaction-diffusion models. We chose to work with a gen- 
eral two- variable reaction-diffusion system which mathe- 
matical description in terms of a nonlinear partial differ- 
ential equation. This equation is written as 



dtV^F2{u,v) + D^\'^v. 



(1) 
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Here u and v represent the concentrations of the reagents; 
Fi(u, v) and F2{u, v) are the reaction functions. The dif- 
fusion coefficients of these rea gent s are Du and Dy. Fol- 
lowing the description in Ref.flj], we define a complex 
ftmction Z = (j)^-\-i4>^, where (j)^ = u—u* and 0^ = v — v* . 
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Here u* and v* are the concentrations of the core of spiral 
wave. 

As pointed out in Ref. [13, the sites of the spi- 
ral waves are just the isolated zero points of the com- 
plex function Z — (f)^ + icf)'^, i.e., the spiral wave field 
in two-dimensional surface, and in general, it is called 
the spiral core. The phase field of spiral wave is de- 
fined as the argument of the complex function Z, i.e., 
0{r,t) — arg(Z) = a.ictan{(j)^ / (j>^) . It is easy to see that 
the zero points of Z are just the phase singular points of 
phase field 6{r, t). By using of topological viewpoints, Liu 
et al.|l6| derived a topological expression of the charge 
density, which is written as 



N 



p{x,y,t) = J2WiS'{r~n{t)], 



(2) 



where Wi is the topological charge of the l-th spiral wave. 
This expression also reveals the topological structure of 
the spiral wave. In our topological theory of spiral wave, 
the topological structure of spiral wave will play an essen- 
tial role. In order to make the background of this paper 
clear, in this section we rewrite the charge density ^ of 
spiral wave as a topological current. 

Now we begin to derive the topological current form 
of the charge density of the spiral core. We know that 
the complex function Z = (f)^ + icj)'^ can be regarded as 
the complex representation of a two-dimensional vector 
field Z = Let us define the unit vector: n° = 

l^(a = 1,2; = = Z*Z). It is easy to see that 
the zeros of Z are just the singularities of n. Using this 
unit vector n, we define a two dimensional topological 
current 

r = J-e'i'^Eabdju'^dkn'', hJ, fc = 0, 1, 2. (3) 

ZTT 

Applying Duan's topological current theory, one can ob- 
tain 



X 



where the Jacobian -D'(^) is defined as 



(4) 



(5) 



The delta function expression ^ of the topological cur- 
rent J* tell us that it does not vanish only when the spiral 
waves exist, i.e.,Z — 0. The sites of the spiral core deter- 
mine the nonzero solutions of J*. The implicit function 
theorem [s^l shows that under the regular condition 



the general solutions of 



a= 1,2, 



(6) 



(7) 



can be expressed as 



x^zi{t), l^l,2,---,N, 



x'^^t. 



From Eq.©, it is easy to prove that 



dx' 
~dt 



(8) 



(9) 



According to the (5-function theory [33| and Duan's topo- 
logical current theory, one can prove that 



N 



dx' 



1=1 



(10) 



in which the positive integer /3; is the Hopf index and 
r]i — sgn{D{^)\^ ) = ±1 is the Brouwer degree. Now 
the density of spiral wave are expressed in terms of the 
complex function Z: 



2tt 

N 



= Pim^"^ - zi). 



(11) 



This is just the charge density of spiral in Ref . [14] . There- 
fore, the total charge of the system given can be rewritten 
as 



„ N N 

1=1 1=1 



(12) 



where Wi is just the winding number of Z around z/, the 
above expression reveals distinctly that the topological 
charge of spiral wave is not only the winding number, but 
also expressed by the Hopf indices and Brouwer degrees. 
For multi-armed spiral wave, the Hopf indices /3; charac- 
terizes the numbers of arms in a spiral wave (/9/ = 1, a 
one-armed spiral wave; = 2, a two-armed spiral wave). 
And the Brouwer degree 77/, whose sign is very impor- 
tant, characterizes the direction of rotation of spiral wave 
(rji = -|-1, the spiral wave rotating clockwise; rji = —1, he 
spiral wave rotating anti-clockwise). It is seen that our 
topological theory provides a more detailed description of 
the inner structures of spiral wave. The topological inner 
structure showed in Eq. (|lip is more essential than that 
in Eq. this is just the advantage of our topological 
description of the spiral wave. 

According to Eq. ([3]) , it is easy to see that the topolog- 
ical current J' is identically conserved, 



a J' = 0. 



(13) 



This equation implies the conservation of the topological 
charge of spiral wave: 



ftp + V • J = 0, 



(14) 
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which is only the topological property of the complex 
function Z. The conservation of the total topological 
charge Q is a very stronger topological constraint on spi- 
ral wave, many important properties of spiral waves due 
to this topological rule. In our following sections, we dis- 
cuss the generating, annihilating at the limit points and 
encountering, splitting, merging at the bifurcation points 
of complex function Z, and it shows that the conserva- 
tion of the total topological charge Q is also valid in these 
processes. 



which is a parabola in — t plane. From Eq. (|20)) we 
can obtain two solutions xXii) and x\(t\ which give two 
branch solutions (world lines of the spiral waves) . If 

We have the branch solutions for t > t* [see Fig. 1(a)]; 
otherwise, we have the branch solutions for t < t* [see 
Fig. 1(b)]. These two cases are related to the origin and 
annihilation of the spiral waves. 



III. THE GENERATION AND ANNIHILATION 
OF SPIRAL WAVES 

As investigated before, the zeros of the complex func- 
tion Z play an important role in describing the topolog- 
ical structures of spiral waves. Now we begin discussing 
the properties of the zero points, in other words, the 
properties of the solutions of Eq. ((T]). As we knew before, 
if the Jacobian 



X 



(15) 



we will have the isolated zeros of the vector field Z. The 
isolated solutions are called regular points. However, 
when the condition (|15p fails, the usual implicit func- 
tion theorem 32] is of no use. The above discussion will 
change in some way and will lead to the branch process. 
We denote one of the zero points as {t* , zl). If the Jaco- 
bian 



^0, 



(16) 



(t%il) 



we can use the Jacobian D^(-) instead of D^(-) for the 
purpose of using the implicit function theorem. Then we 
have a unique solution of Eq. ([7]) in the neighborhood of 
the limit point (t*,z;) 



t^t{x'^), x^^x^{x^), 



(17) 



with t* = t{zl). We call the critical points (t*,i/) the 
limit points. In the present case, we know that 




(a) 




(b) 



FIG 

j2 



1: Projecting the world lines of spiral waves onto 
t) plane, (a) The branch solutions for Ea. (l20|l when 
> 0, i.e., a pair of spiral waves with op- 
posite charges generate at the limit point, i.e., the origin 
of spiral waves, (b) The branch solutions for Ea. l|2Up when 
d^t/ {dx^Y\(^t* ,zi) < 0, i.e., a pair of spiral waves with opposite 
charges annihilate at the limit point. 



~df 



dt 

dx^ 



= 



(18) 



(19) 



Then the Taylor expansion of t = t(x^) at the limit point 
{t*,zi) is 



t-t* 



1 dH 



2 {dx^f 



(20) 



(t*,?l) 



One of the results of Eg. pH)) . that the velocity of 
the spiral waves are infinite when they are annihilating, 
agrees with the fact obtained by Bray [l^ who has a 
scaling argument associated with the point defects final 
annihilation which leases to a large velocity tail. From 
Eg. pS)) . we also obtain a new result that the velocity 
field of spiral waves is infinite when they are generating, 
which is gained only from the topology of the complex 
function Z . 

Since topological current is identically conserved, the 
topological charge of these two generated or annihilated 
spiral pair must be opposite at the limit point, i.e.. 



(22) 
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which shows that = /Si^ and rji-^ = —rji^. One can see 
that the fact the Brouwer degree rj is indefinite at the 
hmit points imphes that it can change discontinuously at 
hmit points along the world lines of the spiral vortices 
(from ±1 to =f1). It is easy to see from Fig. l:when 
> zf,iii^ = ±1; when < zl,rji^ — =Fl. 
For a limit point it is required that D^{^)\^^, 7^ 0. 

As to a bifurcation point [sH], it must satisfy a more 
complex condition. This case will be discussed in the 
following section. 



Case 1(A ^ 0). For A = 4(5^ - AC) > 0, we get 
two different directions of the velocity field of the phase 
singularity for spiral waves 



~dt 



A 



(28) 



which is shown in Fig. 2. It is the intersection of two 
phase singularity for spirals, which means that two phase 
singularity for spirals meet and then depart from each 
other at the bifurcation point. 



IV. BIFURCATION OF THE VELOCITY FIELD 
OF THE PHASE SINGULARITY FOR SPIRAL 
WAVE 



In this section we have the restrictions of Eq. 
the bifurcation points {t*,zi), 



X 



0, 



X 



0, 



jZl at 



(23) 



which leads to an important fact that the function rela- 
tionship between t and x^ is not unique in the neighbor- 
hood of the bifurcation point {t*,zi). It is easy to see 
that 




dt 



(24) 



which under Eq. ipS]) directly shows that the direction of 
the integral curve of Eg. ((24)) is indefinite at (t*,i/), i.e., 
the velocity field of the phase singularity for spirals is 
indefinite at {t*,zi). That is why the very point (i*,i/) 
is called a bifurcation point. 

Assume that the bifurcation point (i*,i/) has been 
found from Eqs.© and (^5]) . We know that, at the bi- 
furcation point (t*,zi), the rank of the Jacobian matrix 
[^] is 1. In addition, according to the (/i— mapping the- 
ory, the Taylor expansion of the solution of <p^ and (jP in 
the neighborhood of the bifurcation point can generally 
be denoted as 

Aix^ - zlf + 2B{x^ - zl){t - t*) + {t- t*f = 0, (25) 



which leads to 



and 



, .dx"^ dx^ ^ 



+ + 



dx^ 



(26) 



(27) 



where A, B, and C are three constants. The solution of 
Eq. ([26)1 or Eq. ([27|) give different directions of the branch 
curves (world lines of the spiral waves) at the bifurcation 
point. There are four kinds of important cases, which will 
show the physical meanings of the bifurcation points. 



FIG. 2: Projecting the world lines of spiral waves onto {x^ —t) 
plane. Two spiral waves meet and then depart at the bifur- 
cation point. 

Case 2{A ^ 0). For A = 4(5^ -AC) = 0, the direction 
of the velocity field of the phase singularity is only one 



dx^ 
~dt 



1,2 



-B 



(29) 



which includes three important situations, (a) One world 
line resolves into two world lines, i.e., one spiral wave 
splits into a spiral pair at the bifurcation point [see 
Fig. 3(a)]. (b) Two world lines merge into one world line, 
i.e., a spiral pair merge into one spiral at the bifurca- 
tion point[see Fig. 3(b)]. (c) Two world lines tangentially 
contact, i.e., a spiral pair tangentially encounter at the 
bifurcation point[see Fig. 3(c)]. 

Case 3(A = 0, C 7^ 0). For A = 4(5^ - AC) = 0, we 
have 



dt 
d^ 



1,2 



B ± - AC 
C 



There are two important cases: (a) Three world lines 
merge into one world line, i.e., three spirals merge into 
a spiral at the bifurcation point [see Fig. 4(a)]. (b) One 
world line resolves into three world lines, i.e., a spi- 
ral splits into three spirals at the bifurcation point [see 
Fig.4(b)]. 
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t' t 

(b) 




t" t 

(c) 

FIG. 3: (a) One spiral wave splits into two spiral wave at the 
bifurcation point, (b) Two spiral waves merge into one spiral 
wave at the bifurcation point, (c) Two world line of spiral 
waves tangentially intersect, i.e., two spiral waves tangentially 
encounter at the bifurcation point. 

Case 4(A = C = 0). Equation(l26l) and give 
respectively 

dx^ (31) 

dt ' dx^ 

This case is obvious, see Fig. 5, and is similar to Case 3. 

The above solution reveal the evolution of the spirals. 
Besides the encountering of the spiral waves, i.e., a spiral 
pair encounter and then depart at the bifurcation point 
along different branch cures [see Fig. 2 and Fig. 3(c)], it 
also includes splitting and merging of spirals. When a 
multi-charged spiral moves through the bifurcation point. 



it may split into several spirals along different branch 
curves [see Fig. 3(a), Fig. 4(b), Fig. 5(b)]. On the contrary, 
several spirals can merge into a spiral at the bifurcation 
point [see Fig. 3(b) and Fig. 4(a)]. 




t" t 

(b) 

FIG. 4; Two important cases of Eq. pO)) . (a) Three spiral 
waves merge into one at the bifurcation point, (b) One spiral 
wave splits into three spiral waves at the bifurcation point. 

The identical conversation of the topological charge 
shows the sum of the topological charge of these final 
spirals must be equal to that of the original spirals at the 
bifurcation point, i.e., 

T.f^hVh^Y.(^hmf (32) 

i f 

for fixed I. Furthermore, from the above studies, we see 
that the generation, annihilation, and bifurcation of spi- 
rals are not gradually changed, but suddenly changed at 
the critical points. 



V. APPLICATION AND DISCUSSION 

Our conclusions can be summarized as follows: First, 
by making use of Duan's topological current theory^ we 
obtain the charge density of spiral waves and the topo- 
logical inner structure of its topological charge, the spi- 
ral waves can be classified not only by their topological 
charges, but especially importantly by their Hopf indices 
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t' t 

(b) 

FIG. 5: Two world lines intersect normally at the bifurcation 
point. This case is similar to FIG. 4. (a) Three spiral waves 
merge into one at the bifurcation point, (b) One spiral wave 
splits into three spiral waves at the bifurcation point. 

and Brouwer degrees. This is more general than usu- 
ally considered and will be helpful as a complement of 
the current description of spiral wave only by the wind- 
ing numbers in topology. Second, the evolution of spiral 
wave is also studied from the topological properties of a 
two-dimensional vector field Z. We find that there exist 
crucial cases of branch processes in the evolution of the 
spiral wave when D{-) ^ 0, i.e., rji is indefinite. This 
means that the spiral waves are generate or annihilate 
at the limit points and encounter, split, or merge at the 
bifurcation points of the two-dimensional vector field Z, 
which shows that the spiral wave system is unstable at 
these branch points. Third, we found the result that 



the velocity of spiral core is infinite when they are anni- 
hilating or generating, which is obtained only from the 
topological properties of the two-dimensional vector field 
Z. Forth, we must pointed out that there exist two re- 
strictions of the evolution of spiral waves. One restriction 
is the conservation of the topological charge of the spiral 
waves during the branch process [see Eqs. and ([5^ ]. 
the other restriction is that the number of different di- 
rections of the world lines of spiral waves is at most 4 at 
the bifurcation points [see Eqs. ([26]) and ([27]) ]. The first 
restriction is already known, but the second is pointed 
out here for the first time to our knowledge. We hope 
that it can be verified in the future. Finally, we would 
like to point out that all the results in this paper have 
been obtained only from the viewpoint of topology with- 
out using any particular models or hypothesis and do 
not require the knowledge of the actual solution of the 
reaction-diffusion equation. 

In this paper, we give a rigorous and general topo- 
logical investigation of spiral wave. This work can be 
applied in several ways. First, topology and geometry of 
two-dimensional manifold in which the spiral lived can 
affect the dynamics of spiral wavefl, Q, therefore, it is 
very interesting to extend our work to a curved surface 
which have complex topology and geometry, and study 
the interaction between the topology of spiral wave and 
the topology of curved surface. Second, our theory can be 
easily extend to the three dimensional scroll wave, which 
is just a three dimensional analog of spiral wave, in this 
situation, a more complex topology related to knots must 
be considered d^. Third, our work provide a more de- 
tailed description of spiral wave, the detailed information 
of spiral can applied to the study of dynamics of spiral 
wave. These issues are very worthwhile to consider and 
will be investigated in our further works. 
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